Coulomb Blockade due to Quantum Phase-Slips Illustrated with Devices 
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In order to illustrate the emergence of Coulomb blockade from coherent quantum phase-slip 
processes in thin superconducting wires, we propose and theoretically investigate two elementary 
setups, or "devices". The setups are derived from Cooper-pair box and Cooper-pair transistor, so 
we refer to them as QPS-box and QPS-transistor, respectively. 

We demonstrate that the devices exhibit sensitivity to a charge induced by a gate electrode, this 
being the main signature of Coulomb blockade. Experimental realization of these devices will unam- 
biguously prove the Coulomb blockade as an effect of coherence of phase-slip processes. We analyze 
the emergence of discrete charging in the limit of strong phase-slips. We have found and investigated 
six distinct regimes that are realized depending on the relation between three characteristic energy 
scales: inductive and charging energy, and phase-slip amplitude. For completeness, we include a 
brief discussion of dual Josephson-junction devices. 



I. INTRODUCTION 

The interplay between superconductivity and Coulomb 
interactions in mesoscopic superconducting circuits is a 
subject of active research for more than 20 years 1 - - — . A 
generic device that has been widely studied and used is 
the Cooper-pair box (CPB) 3 ' 6 ' 7 . It consists of a super- 
conducting island to store discrete charges that is con- 
nected to a bulk electrode by means of a tunnel Joseph- 
son junction. This enables coherent transfer of Cooper 
pairs between the island and electrode. A gate electrode 
capacitively coupled to the island induces a continuous 
charge q. The energies of quantum states of the CPB are 
periodic functions of this induced charge. This charge 
sensitivity reveals the charging states of the device, those 
with a well defined number of excess Cooper pairs, and 
opens up the possibility to create and control their quan- 
tum superpositions. For instance, the CPB can be oper- 
ated near the point where two charging states are approx- 
imately degenerate forming a qubit basis. The Joseph- 
son tunneling lifts the degeneracy between these charging 
states and enables their superpositions. 8 

CPB and similar devices based on Coulomb blockade, 
like flux qubits^, have been used to realize viable qubit 
schemes. Among other achievements, coherent control of 
quantum states 1 ^, Rabi oscillation a 11 ' 12 , successful DC 
readou t 13 i 14 and RF coupling of multiple junctions have 
been demonstrated. 

The studies of phase-slip processes 1 ^ in superconduct- 
ing wires have a long and spectacular history. During a 
phase-slip process, the superconducting order parameter 
passes zero at a certain moment of time and at a certain 
position in the wire. The phase difference between the 
wire ends changes by 27r. In accordance with Josephson 
relation, this gives rise to a voltage pulse across the wire. 
It has been established more than 40 years ago 1 ^ that the 
residual resistance of the wires below the superconduct- 
ing transition is due to thermally activated phase-slips»il 

At sufficiently low temperatures, the quantum fluctua- 
tions should supersede the thermally-activated ones. The 



manifestations of these quantum phase-slips in ultra-thin 
resistive wires have been actively investigated for the last 
10 yearsi^r— . It has been suggested that the switch- 
ing from superconducting to normal state in current- 
driven wires is caused by individual quantum phase-slips 
event o 21 ' 22 . These experiments implemented resistive 
measurements where the quantum coherence between in- 
dividual phase-slip events seems to be destroyed by ac- 
companying dissipation. If such coherence is preserved, 
the manifestations of phase-slips is qualitatively differ- 
ent: they change the characteristics of the ground and 
excited quantum states^ It has been proposed 2 ^ that 
the coherent phase-slips should be observed and studied 
in non-driven devices of qubit type. Recently, the ef- 
fect of coherent quantum phase-slips has been observed 
in Josephson chains that are in many respects similar to 
superconducting wires 2 ^. 
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FIG. 1. a. Cooper-pair box(CPB): generic Coulomb-blockade 
superconducting device, b. In quantum phase-slip box (QPS- 
box), the tunnel junction is replaced with a superconducting 
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FIG. 2. Energy spectrum of CPB versus induced charge q 
in the limit of large charging energy. The discrete charging 
states give rise to a standard 2e-periodic pattern of crossing 
parabolas. The spectrum of phase-slip devices considered in 
the article follows the same pattern in the limit of large phase- 
slip amplitudes. 



It is clear that sufficiently well-developed coherent 
phase-slips should eventually lead to Coulomb blockade 
in the wire. This statement is fascinatingly counterintu- 
itive: the superconducting wire is one of the best con- 
ductors possible, while Coulomb blockade requires isola- 
tion. However, the necessity of Coulomb blockade and 
accompanying isolation in this regime directly follows 
from basic arguments that involve duality between the 
charge and phased. In a nutshell, the argumentation 
is as follows. The coherence of Cooper-pair tunneling 
events leads to a zero-voltage states at currents below 
a critical value. Interchanging charge and phase, current 
and voltage, we prove that the coherence of the phase-slip 
events should lead to a zero-current state below a critical 
voltage. This is essentially Coulomb blockade and isola- 
tion. Indeed, there are experimental data that indicate 
a possible crossover to insulating behavior in ultra-thin 
wiresj 19 ' 21 However, the exact conditions and mechanism 
of this crossover are still subject to debated 

The main motivation of this article is to facilitate 
an unambiguous experimental proof of Coulomb block- 
ade due to phase-slips. In this context, we note that 
Coulomb blockade is usually accompanied by periodic 
charge sensitivity. The observation of charge sensitivity 
in setups where the tunnel barriers are replaced by un- 
interrupted superconducting wires would constitute the 
proof required. The experimental attention to charge 
sensitivity is presently insufficient, though a very recent 
communication 28 reports indirect observation of charge 
sensitivity by its effect on the decoherence in fluxonium 29 
qubit. We have also recently learned of unpublished re- 
sults of the authors of Reb 2 ^ that demonstrate the gate- 
voltage effect. In this article, we analyze the problem 
at elementary level by introducing simple phase-slip sys- 
tems (we call them devices) where this charge sensitivity 
can be observed, and discuss the conditions for this to 
happen. 



The devices are obtained from the generic Cooper-pair 
box (CPB) and Cooper-pair box with two leads (Cooper- 
pair transistor, CPT) by replacing the tunnel junctions 
with the superconducting wires subject to phase-slips. 
We ascribe all the capacitance to the island, thus dis- 
regarding geometric capacitance of the wires. We also 
disregard quasiparticle effects assuming that the super- 
conducting gap by far exceeds the temperature. These 
assumptions allows us to account for the phase-slips in 
the framework of a zero-dimensional phenomenological 
model of phase-slip junction^ (we discuss the relation 
of this model and microscopic theory of phase-slips in 
the Appendix). We call the devices QPS-box and QPS- 
transistor. 

Let us note that for the devices proposed, the Coulomb 
blockade and associated charging states do not occur in 
the wire. It would be wrong to assume that the wire is 
broken into a chain of weakly connected Coulomb islands. 
Indeed, the assumption of vanishing self-capacitance in 
fact forbids any charge accumulation in the wire. Rather, 
the collective state developed in the wire as a whole pro- 
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FIG. 3. a. Cooper-pair transistor (CPT). The dependence 
of its energy on external phase $ (flux sensitivity) gives rise 
to the supercurrent that can be modulated by changing q = 
CgVg. b. Quantum phase-slip transistor. 
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vides the isolation between the wire ends. Qualitatively, 
the wire as a whole acts as a single tunnel barrier. 

We study the energy levels in these devices, with em- 
phasis on the ground state. The operational limit of the 
devices is set by the value of the effective impedance 
7 = (El/Ec) 1 ^ 4 - the ratio of the inductive and capac- 
itive energy. In the limit of large impedance, the phase 
fluctuations are dominant, while for small impedance 
regime the charge fluctuations become relevant. The val- 
ues of L and C are set by fabrication and they have a 
wide range of possible values. In order to make quali- 
tative and quantitative predictions about the devices we 
compare the corresponding inductive and charging ener- 
gies with the phase-slip amplitude. 

We demonstrate that these energies exhibit charge sen- 
sitivity (for the boxes) and combined flux/charge sen- 
sitivity (for the transistors). Moreover, in the limit of 
large phase-slip amplitude, the phase-slip devices can be 
directly mapped on their tunnel-junction counterparts, 
and exhibit the standard pattern of " crossing parabolas" : 
charge-sensitive discrete charging states. 

For completeness, let us make a remark about the ef- 
fect of the random offset charges on our devices. With- 
out phase-slips the proposed devices are just linear and 
not affected by any charge. However, the charge sen- 
sitivity brought about by the phase-slips implies that 
the devices are affected not only by the gate voltage but 
also by the random offset charges most likely present in 
the substrate. In general, it is a task of the fabrication 
technology to minimize the effect of offset charges. We 
expect, that the effect of the offset charges on the de- 
vices under consideration would be the same as that on 
any other Coulomb blockade system like SET transistor, 
quantum dot, superconducting qubit etc. At sufficiently 
low substrate temperature the random offset charges re- 
main the same and can be compensated for by a shift in 
the gate voltage. If the devices are realized with metal 
wires (rather than with Josephson junction chains) there 
is a chance that the offset charges are efficiently screened 
by the metal leads. Therefore, such realization may be 
better than that involving tunnel junctions where the 
offset charges are present in the insulating layers of the 
junctions. 

The paper is organized as follows. Section ITT1 describes 
in detail the devices under consideration and establishes 
their Hamiltonians. Section Mil is devoted to QPS-box. 
We consider first the phase-slips perturbatively, then an- 
alyze the crossover to Coulomb blockade in the regimes 
of large, small and intermediate impedance. Section IIVI 
details the QPS-transistor with emphasis on flux sensi- 
tivity specific for this device. In Section|V]we discuss the 
Josephson-based devices that are dual to QPS-box and 
QPS-transistor. We conclude in Section IVT1 



II. DESCRIPTION OF THE DEVICES 

Let us recall two generic devices that exemplify the 
manifestation of Coulomb blockade in superconducting 
circuits. They are made by connecting a superconduct- 
ing island with either one or two superconducting leads. 
An isolated island supports discrete charges. The im- 
portant part of the setup is the gate electrode that is 
not electrically connected to the island but, by means 
of capacitive coupling, induces charge q on the island. 
We will refer to these two devices as to Cooper-pair box 
(CPB, Fig. Ha.) and Cooper-pair transistor (CPT). The 
latter term is less conventional: we use it because the su- 
percurrent through the device does depend on the gate 
voltage, this being a transistor effect. Besides, the setup 
reminds much that of normal-metal Single-Electron Tun- 
neling Transistor (SET)22. 

Our idea of introducing the phase-slip devices is to 
replace the tunnel junctions in above setups with thin 
superconducting wires. In this way, we come to the se- 
tups of QPS-box and QPS-transistor. Healthy reasoning 
is that the wires would short-circuit superconducting is- 
land to the lead or leads: unlike the tunnel junctions, the 
wires are not expected to provide isolation required for 
Coulomb blockade phenomena. In this case, the charge 
induced to the island should have no physical effect. This 
is indeed true if coherent phase-slips in the wires are dis- 
regarded. In this case, both devices are just linear electric 
circuits. 

It is the main goal of the present manuscript to show 
in detail that the coherent phase-slips induce charge sen- 
sitivity in both devices. Moreover, for sufficiently large 
phase-slip amplitudes, the phase-slip devices are eventu- 
ally identical to CPB and CPT, respectively, and exhibit 
almost pure charging states like shown in Fig. [3] For 
QPS-transistor, the charge sensitivity is combined with 
the flux sensitivity, so one can observe a transistor effect. 

All devices are simple to describe quantum- 
mechanically: this involves a Hamiltonian with several 
degrees of freedom only. In the rest of this Section, 
we give these Hamiltonians as the basis for further 
consideration. 



A. CPB versus QPS-box 

The Hamiltonian of the CPB consists of charging and 
Josephson terms, 

^CPB = E c(Q-^) 2 -^JCOS(^), (I) 

o 2 

where Eq = (c+c ) lii ^ ne charging energy involving the 
total capacitance of the island. The Josephson energy 
can be expressed in terms of the conductance of the tun- 
nel junction G, Ej = -j , where Gq is the conductance 
quantum, and A the superconducting gap. The operator 
Q is the charge stored on the island measured in units of 
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Cooper-pair charge 2e, <f> is the operator of the supercon- 
ducting phase difference between the island and the lead. 
These two operators satisfy the canonical commutation 
relations [Q,4>] — —i- 

This Hamiltonian is undefined unless we specify the 
space in which the variables Q,(f> are defined. For CPB, 
the phase is defined in the interval (— it, n) so that the 
wave function is periodic in <fr, ip(<fi) = ip((f> + 27r). This 
assumption makes the charge variable discrete, Q = N, 
N being integer numbers spanning a countable set of 
states \N). The Hamiltonian in charge representation 
then reads 



H CPB = Y, e c(n- \N) (N\ 

N 



Ej 



(\N) (N + l\+h.c). 



(2) 



The QPS-box includes a superconducting wire. If we 
restrict ourselves to the energies lower than the inverse 
time of propagation of electric excitations along the wire, 
we can neglect the details of spatial distribution of the 
order parameter in the wire and characterize its quantum 
state by a single quantum variable <f> that correspond to 
the phase drop across the wire. Neglecting the details 
of the spatial distribution implies that we the disregard 
geometric capacitance of the wire. More precisely, we 
ascribe all the capacitance to the island. We stress that 
in distinction from CPB the phase variable is extended 
being defined in the interval (—00,00). One can regard 
this variable as consisting from a compact phase of the 
island and discrete number of phase windings p along the 
superconducting wire. Correspondingly, the charge is a 
continuous variable, and Coulomb blockade is not evi- 
dent. Without phase-slips, the wire can be regarded as a 
linear inductor of inductance L. Together with the island 
capacitance, this gives an LC oscillator. The phase-slips 
change the winding number p by ±1. Since we neglect the 
details of spatial distribution, it does not matter where in 
the wire a phase-slip would occur. The tunneling between 
different p can be thus described by a single amplitude. 26 
We present more microscopic details in the Appendix. 

These assumptions define the Hamiltonian of the QPS- 
box. It consists of a part that describes a linear circuit, 
in our case, an LC-oscillator, and phase-slips 
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(q 



2e 
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4 ' 



Ha 



(3) 



phase variable by ±2it: 

H s V#) - -E s V# + 2tt) - E s - 2tt), (4) 

E$ being the phase-slip amplitude. In the charge repre- 
sentation this term takes a simple form of a cos-potential: 



H s = -2E s cos(2irQ). 



(5) 



The charge sensitivity is entirely due to the phase-slip 
term: one can make this explicit by shifting the charge 
variable by q/2e. The induced charge q disappears from 
the oscillator term while the phase-slip amplitudes ac- 
quire phase factors: 



H s tp{4>) = -E s e~" q/e ip(4> + 2tt) - E s e™ q/e ip(4> - 2tt). 

(6) 

One can say that the induced charge affects the in- 
terference of phase-slips of two opposite directions. It 
is a matter of choice whether Eq. [4] of Eq. [6] defines 
the phase-slip operator since the difference is the shift 
in charge space. Depending on the choice, the charging 
energy reads either Ec(Q — q/2e) or EcQ- 

We notice that the QPS-box without phase-slips (Es — 
0) is a LC oscillator with the frequency 



h^E L E c 



and the spectrum 



E n = hujQ n + 



1 



n being number of oscillator quanta. The oscillator is 
characterized by an important parameter 7 defined as 



7 : 



that measures the effective impedance Z of the oscilla- 
tor in quantum units. This parameter can be readily 
expressed through the ratio of the charging and induc- 
tive energies: 7 2 = j^. In terms of the creation- 
annihilation operators b, of the oscillator the phase-slip 




where 



Hamiltonian reads Hs — —Es(e 



y(fct- b ) 



h.c). 



Ec — 



c 



Er = 



$0 

2tt 



are the charging energy corresponding to charging the 
capacitor C and respectively the inductive energy corre- 
sponding to the inductor L, $0 being the flux quantum. 

The action of the phase-slip Hamiltonian adds ±1 to 
the winding number, that is, shifts the wavefunction in 



B. CPT versus QPS-transistor 

Physically, the CPT is a device that is more functional 
than the CPB. It can conduct the superconducting cur- 
rent that is affected by the induced charge, this presents 
a convenient way to detect the charge sensitivity. How- 
ever, at Hamiltonian level the devices are described by 
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the same Eq. [3J The only difference is that the Cooper- 
pair tunneling to/from the island can proceed through 
two junctions, both adding to the tunneling amplitude, 



Er = E. 



7.1 



Ej.ie 



t'i> 



The superconducting phase difference between the leads 
$ affects the interference of the two amplitudes. We note 
that a practical way to realize such phase bias is to em- 
bed the device into a superconducting loop penetrated 
by magnetic flux so that we interchangeably refer to this 
parameter as to phase or flux. For instance, we talk 
about flux sensitivity while referring to the dependence 
of energy levels on $. 

In contrast to this, the Hamiltonian description of the 
QPS-transistor (Fig. Ob) is more complex than that of 
the QPS-box. The point is that there are two wires in 
the device, each with its own winding number. This gives 
rise to two extended phases </>i.2- However, there is a con- 
strain on these phases. They should sum up to the overall 
phase drop over the device that is fixed externally. It is 
convenient to restrict the external phase to the interval 
(— 7r,7r). In this case, the constrain reads 



b 2 = $ + 2irp, 



(7) 



integer p being the total winding number in both wires. 

Let us introduce two continuous charge variables Qi,2 
that are canonically conjugated to these phases. With 
this, the Hamiltonian of the QPS-transistor reads 

El,1 22 



-Hpst = 



E L,2 22 



E C {Q 2 - Qi - 



(8) 



+ Hs,i + Hs,2, 

where the phase-slip Hamiltonians are 

H S,\(2) = -2£ , s,i(2)C0s(27rQ 1( 2)). 

We recognize the inductive terms, proportional to the 
squares of the phase drops </>i,02 across each inductor. 
The total charge accumulated in the island is the differ- 
ence of the charges Q2,Qi passing the wires, this gives 
the form of the charging energy. The last two terms 
present the phase-slips that shift the corresponding phase 
drops by ±2tt. The Hamiltonian can be written in several 
equivalent representations. In charge representation, the 
constrain (0 implies the the following periodicity condi- 
tion on the wave function, 

#2i,Q 2 ) = ^(Qi + i,<92 + i) 

while the Hamiltonian reads 



Hpst = 



E 



. d 



E 



L,2 



. d 
'dQ~2 



E c Q 2 -Qi-^- 
2E s ,i cos(2ttQi) - 2E S ,2 cos(2irQ 2 



(9) 



However, the properties of the QPS-transistor at 
Es,i,Es,2 are easier to understand in phase representa- 
tion. In this case, the wavefunction is defined at a series 
of lines numbered by p and parametrized with contin- 
uous <f>\. At each line, the wavefunctions are that of a 
harmonic oscillator with the frequency 



f^o = \JE c (E Lil + E Lt2 ). 

The ground state energies of these oscillators are different 
on different lines. They correspond to minimum induc- 
tive energy at a given total winding number p, 



E„ 



El,i ■ El : 



p 4(E L ,i + E L , 2 ) 
Therefore the spectrum reads 



(2irp + $) 2 



E n , p = haj(n+^j + E p . 



(10) 



It is important to recognize that the oscillators at differ- 
ent p arc shifted with respect to each other, their equi- 
librium positions being given by 



4%) 



El, 



El,i + El ,2 



(2irp + $) 



In this representation, the phase-slip Hamiltonians ei- 
ther increase or decrease the total winding number, while 
Hs,i also shifts 4>i by ±27r, 



H s ,ii>{4n,p) = -2E s ,iV(<ln +2tt, P +1) 

-2E$ tl *(<l> 1 -2*,p-l), (11) 

Hs^i^p) = -2E s ^ uP + 1) - 2E* S 2 ip(<pi,p- 1) 

(12) 

We define the parameter 7 that measures the effective 
impedance of the oscillator j 2 = yj (El,i + El.2)/Ec- 

In contrast to the QPS-box, the QPS-transistor ex- 
hibits both flux and charge sensitivity that makes it po- 
tentially useful for measurements. 

III. QPS-BOX 

In this Section, we study the ground state and low- 
energy states of QPS-Cooper pair box, with emphasis 
on the charge sensitivity of their energies. The charge 
sensitivity appears already in the limit of E$ — > 
as a first-order perturbation correction to the oscilla- 
tor levels. Upon increasing Eg, the charge sensitivity 
increases and eventually the low-lying states follow the 
standard Coulomb-blockade pattern. The crossover to 
the Coulomb blockade follows different scenarios depend- 
ing on the effective impedance of the oscillator I/7 2 . 
We discuss the limits of large and small impedance and 
present numerical illustrations for these limits as well as 
for the case of intermediate impedance. 
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FIG. 4. The magnitude of the charge-sensitive first-order cor- 
rection (Eq. I13[) to the first four energy levels n = 0, ... 4, 
versus 7 2 . 

A. First-order corrections 

In the limit of E$ <C Ec , El it should be possible 
to treat the phase-slip term H$ as a perturbation. The 
unperturbed system at E$ = is nothing but an LC 
oscillator, and the first order correction is given by the 
diagonal matrix element of H$, E^ = (n\Hs\n), with 
respect to these states. The correction gives the charge- 
sensitive part of the energy and reads 

E$p = -2.B5COS (— ) exp(-7rV/2) iFi[-n, 1,7V 2 ]. 

(13) 

Here \F\ stands for the confluent hypergeometric func- 
tion of the first kind. 

This expression is valid at any value of the effective 
impedance. In Fig. 2] we plot this correction for the first 
four eigenenergies versus 7 2 . For the ground state n = 0, 
the correction reads 

^=-2^cos(^)ex P (4^. (14) 

we we use 7 2 = yj Jp^. 

The first-order corrections have a characteristic 
cos(7rg/e) dependence on the induced charge. Despite 
their relatively small magnitude, they can be revealed 
by driving the oscillator with a.c. gate voltage at a 
frequency close to uioM- The first-order corrections are 
exponentially suppressed in the limit of small effective 
impedance, 7 > 1. To see that, let us estimate a typ- 
ical spread of the wave function of the ground state in 
phase space, that is, quantum fluctuation of phase. For 
the case of small effective impedance, this fluctuation is 
small, ~ I/7. The first-order correction due to phase- 
slips is given by the overlap of this ground state wave 
function with its copy shifted by 2n. The small spread 
leads to exponentially small overlap and thus to the ex- 
ponentially small correction. 

Comparing these first-order corrections with oscillator 
energy differences Hluq suggests that the crossover to well- 
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FIG. 5. QPS-box: Small impedance regime - Bloch states: 
The energies versus quasi-phase x- a: Strong PS:Es = 
6.25-Ei. b: Moderate PS: E s = 2.5E L . c: Weak PS: 
E S = 0.25E L . 

developed Coulomb blockade takes place at Es — hu>o for 
7 < 1 and at exponentially large Es for 7> 1. This esti- 
mation is too simplistic: we will consider below the cases 
of small and large effective impedance to show where and 
how the crossover actually happens. 

B. Small impedance regime 

In this case, El ^> Ec and thus 7> 1. To understand 
the specifics of the regime, let us first set the charging 
energy Ec to zero. The Hamiltonian in charge represen- 
tation reads 

is translational invariant in Q-space, and is equivalent 
to a Hamiltonian of a quantum particle a periodic po- 
tential. Owing to translational invariance, the eigenen- 
ergies do not depend on q showing no charge sensitiv- 
ity. The eigenfunctions are Bloch states labeled by x> 
i/j(Q + 1) = e tx i/j(Q). We will call the parameter \ the 
quasiphase, since it is similar to the variable (f> and plays 
the role of quasi-momentum for common Bloch states in 
periodic solids. 

The properties of Bloch states are determined by com- 
petition of the energy scales Es and El- If El 3> Eg, 
they form nearly parabolic subbands (see Fig. [S]c). The 
energy correction to the ground state appears in the sec- 

(2) E 2 

ond order in Eg, Eq oc -jf-. ) Note that in distinction 
from the first-order correction, Eq. [14J this is not expo- 
nentially suppressed. 

In the opposite limit of El <C Es, we have a well- 
developed periodic potential in Q-space. The lowest 
Bloch subbands correspond to quantized energy levels 
in equivalent potential wells. Near the potential mini- 
mum, the energies of these levels are those of an effec- 
tive oscillator, E n = n^/E~s~E~L(2n +1). Their dispersion 
comes about the tunneling through potential barriers of 
the height ~ Es and is exponentially suppressed(see Fig. 
[5] a). For the lowest subband, where the suppression is 
the strongest, 

E = -\E 8 \ + -k^/ElEs + 2A cos(x), 

A = 8{E 3 S E L ) 1/4 ■ exp(-%J§p\ being the ampli- 
tude of tunneling between the lowest-energy states in the 
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FIG. 6. The renormalized inductive energy E L versus the 
phase-slip amplitude Es- It becomes exponentially sup- 
pressed at Es S> El- 



neighboring wells. 

Let us now resume our analysis of QPS-box and con- 
sider small but finite charging energy Ec- Thereby we 
lift the degeneracy of the potential minima that would 
in principle lead to charge sensitivity. Let us concentrate 
on the lowest energy eigenstates. Those should originate 
from the Bloch states with lowest energy. Their disper- 
sion can be approximated as const + E L /A\ 2 , where we 
expand in the quasiphase \ and introduce an effective 
inductive energy, E* L . The latter is given by the second 
derivative of the lowest energy band energy with respect 
to x an< i represents the renormalization of inductance by 
the phase-slip processes. We plot this quantity versus 
E S /E L in Fig. m 

With this, the lowest energy states can be approxi- 
mated by an effective Hamiltonian in quasi-phase repre- 
sentation, 



4 A 



E, 



c 



d\ 2 



(15) 



which is one of a harmonic oscillator with the renormal- 
ized oscillation frequency ^JW^Ecjh. These states are 
not charge-sensitive in this approximation of a renormal- 
ized oscillator. We expect the charge sensitivity to set on 
only if the approximation breaks down. 

The validity of the approximation can be estimated 
by comparing the spread of the ground state wave func- 
tion in quasi-phase space (~ (Ec/EI) 1 / 4 ) and 2n, the 
scale at which we expect the dispersion of Bloch states 
to deviate from the quadratic law. We thus expect the 
approximation to break down, and the charge sensitiv- 
ity to set on, when Ec — E* L . Since in the limit under 
consideration El ^> Ec, this can only become possible if 
the renormalized inductance is strongly suppressed. This 
requires Es >• El and yields E* L m A . Comparing A 
and Ec, we obtain the estimation for Es at which the 
ground state becomes charge-sensitive, 



E* s 




In this limit, we can restrict our consideration to the 
states of the lowest Bloch subband assuming that the 
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FIG. 7. Small impedance regime (7 = 1.91): the excitation 
energies versus Es/El for two values of the induced charge 
of a: q/e = and b: q/e = 1. The "sticking together" 
of the excitation energies at sufficiently large Es indicates 
emergence of the charging states. 




FIG. 8. The charge sensitivity in the small impedance regime 
(7 = 1.91) for a: E S /E L = 5.48, b: E S /E L = 8.22 and c: 
E S /E L = 10.96. 



energy is small in comparison with the subband splitting, 
E <C \JE L E S > E C - This is equivalent to taking into 
account only the lowest-energy state in each potential 
well. We number these states with N and arrive at the 
effective Hamiltonian 



H = -AoY,(W)(N + l\ + \N- 



N 



l)(N\)+E c (N-q/2e) 2 
(16) 

which appears to be the same as the Hamiltonian Eq. [5] 
of the Cooper-pair box with Ej replaced with 2Ao. The 
energy levels follow the standard Coulomb blockade pat- 
tern provided Ec ^> Aq. Therefore, we have proven the 
equivalence of CPB and QPS-box in the limit of large Es 
and the emergence of well-developed Coulomb blockade 
in the QPS-box. 

We complement this analytical consideration valid in 
the limit 7 3> 1 by numerical calculations at finite value 
7 = 1.91. In Fig. [7] we plot the energies of several excited 
states counted from the ground state versus Es / El ■ To 
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illustrate the charge sensitivity, the plots are made at 
q/e = and q/e = 1. We see how the states evolve from 
equally separated oscillator levels at E$ — to charging 
states at E$ that follow a typical Coulomb blockade de- 
generacy pattern (at q/e = only the excited states are 
doubly degenerate while at q/e = 1 the ground state is 
doubly degenerate as well). We see that the higher ex- 
cited states become more charge-like at smaller values of 
Es, due to quadratic dependence of the charging energy 
on the state number. This is also illustrated in Fig. \E\ 
where we plot the energies of the ground and first two 
excited states versus induced charge upon the increase of 
E s - 
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FIG. 9. Excitation energies in the large impedance regime 
(7 = 0.32) for a: q/e = and b: q/e = 1. 



If Es <C Ec, the energy has a single minimum, and 
exhibits weak charge sensitivity E(q) = —2Es cos(irq/e), 
in agreement with Eq. Q31 At q/e = 1 and criti- 
cal value of Es — Ec /8tt 2 this minimum splits into 
two. More minima emerge upon increasing Es- Fi- 
nally, at Es ^> Ec, a large number (~ irEs/Ec) of low- 
energy minima are pinned by the oscillatory potential to 
Qm — N. Their energy is contributed by the charging 
energy term only and reproduces the Coulomb blockade 
pattern En = Ec(N—q/2e) 2 . We thus conclude that the 
crossover to the well-developed Coulomb blockade takes 
place at Es — Ec ■ The crossings of the energy levels are 
not avoided in this approximation. 

Let us take into account the final value of El- This 
leads to quantization of energy levels around each mini- 
mum that is described by an effective oscillator Hamilto- 
nian, 



E c {Qm — QY 



El_cP_ 
4 <9Q 2 



where E c = 2E"(Q m ) is the charging energy renormal- 
ized by the phase-slips. Since El is small, the the dis- 
tance between the levels is smaller than the scale of the 
potential E(Q). At Es >• Ec we deal with almost iden- 
tical potential wells of depth AEs- The oscillatory mode 
in each well gives rise to a system of equidistant levels 
separated by 2-k\JEsEl- Finite El also induces tun- 
neling between the neighboring wells that removes the 
degeneracy at the level crossings. 

We distinguish two cases in the regime of well- 
developed Coulomb blockade Es 3> Ec- At sufficiently 
small Es, the effective oscillator frequency 2-k\[E~s~El is 
much smaller than the Coulomb energy. In this case, the 
lowest-energy states are charging states (labeled by N) 
with n oscillator quanta. If one neglects the exponen- 
tially small tunneling between the wells, the energies of 
these states are given by 



C. Large impedance regime 

In this limit, Ec >• El and 7 < 1. To start with, 
let us disregard the inductive energy. The Hamiltonian 
is diagonal in charge representation and reads 



H = E ( 



(q 



2e 



2E s cos(2irQ). 



(17) 



The stable states are associated with the minima of 
E(Q). The positions of minima Q m are determined from 
the equation 



E c _ sin(27rQ m ) 

— — = — Z7T — — g — 

Es Qm — 2e 

the corresponding energies being given by 



(18) 



E m — Ec ( {Qm — ^-) 2 + — 

2e 7T 



_ Q_ 

m 2e 



cot(27rQ) (19) 



E(N, n) = Tr V / E s E L (2n + 1) + E C {N - q/2ef 

Tunneling between the wells leads to a fascinating picture 
of avoided crossings between the energy levels that differ 
in N and n (see Fig. (TUJc) 

In the opposite case Es ^> E c /El, the first excited 
state of the oscillator lies much higher than many charg- 
ing states. The tunneling in this case is between the 
ground states in each potential well. The low-energy 
states are described by the same Hamiltonian Eq. \W\ 
as in the regime of small impedance. However, in dis- 
tinction from the small impedance regime, Ao is always 
parametrically smaller than Ec- 

The above analytical analysis is confirmed by numeri- 
cal calculations at 7 = 0.32. In Fig. [9] we see the evolu- 
tion of the excitation energies from those of the original 
oscillator to those of charging states (seen as horizontal 
lines at sufficiently large Es) or charging states with ex- 
tra number of quanta of the renormalized oscillator (seen 
as rising lines owing to Es dependence of the oscillator 
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FIG. 10. The charge sensitivity in the large impedance regime (7 = 0.32) for a: Es/Ec = 0.02, b: Es/Ec = 0.1 and c: 
Es/Ec = 0.3. 




FIG. 11. The evolution of E(Q) at q/e = 1 upon increasing 
the phase-slip amplitude for a: Es/Ec = 0, b: Es/Ec = 0.5 
and c: Es/Ec = 10. The multiple potential minima formed 
correspond to discrete charging states. 
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FIG. 12. Excitation energies in the intermediate impedance 
regime (7 — 1) for (a)g/e = and (b) q/e = 1. 
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FIG. 13. The charge sensitivity at the intermediate 
impedance (7 = 1) for a: Es/fiujo = 1-0, b: Es/ftu>o = 2.0 
and c: Es/hujo = 4.0. 



frequency). The charge sensitivity is illustrated in Fig. 
ITOl . We plot a number of low-lying states versus q at in- 
creasing values of Es- We see the formation of charging 



and oscillator states and progressive reduction of anti- 
crossings between those. 

D. Intermediate impedance regime and summary 

To investigate the regime of intermediate impedance, 
we have numerically computed the cigcncnergies of the 
Hamiltonian Eq. [3] at 7 = 1. The results for excitation 
energies are presented in Fig. [T3] 

As in previous regimes, we observe a crossover from 
equidistant levels of the original oscillator to the charging 
states. The natural scale of for this crossover is neither 
El nor Ec, rather, this is the frequency of the origi- 
nal oscillator hu>o = y/E^Ec- As in the regime of large 
impedance, the charging states are augmented with quan- 
tized oscillations around the minimum of each potential 
well: excitations with n ^ are manifested as the curves 
rising with increasing E$- In distinction from the large 
impedance case, the energy of these excitations slightly 
exceeds the energy of the lowest charging states once they 
are formed. We note that for small impedance the ex- 
citation energies typically decrease with increasing E$ 
while for large impedance they typically grow. In in- 
termediate case, we see that many excitations just keep 
approximately the same energy while converted from an 
oscillator level to a charging state. 

Charge sensitivity (Figll3p shows a pattern similar to 
that in small impedance regime. This proves that at 
sufficiently large E$ in all cases we reach a well-developed 
Coulomb blockade. 

We summarize the regimes of QPS-box in Fig. [TJJ 
This is a log-scale diagram with Es on the vertical axis 
and 7 2 on the horizontal one. Two crossing thick lines in- 
dicate El — Es and Eq — E$, correspondingly. Below 
the lines, Es <C max(i?£, Ec) phase-slips can be con- 
sidered perturbatively (region I). At the lines, the main 
effect of phase-slips is the renormalization of either capac- 
itance (region II) or inductance (region III) of the original 
oscillator. Right above the line Es — Ec we have a re- 
gion IV were the charging states are developed and each 
is accompanied by a set of closely-spaced oscillator levels. 
At the dashed line Es — E c /El, the effective oscillator 
frequency becomes comparable with the charging energy. 
In the region V, the QPS-box is described by the CPB 
Hamiltonian Eq. [TBI At the dashed line, Ec — A . Fi- 
nally, the region VI corresponds to the standard pattern 
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FIG. 14. Distinct regimes in QPS-box. The lines in his 
schematic log-log plot present parameter regions in the space 
of Es and effective impedance I/7 2 . I. Perturbative regime. 
II. Oscillator with renormalized capacitance. III. Oscillator 
with renormalized inductance. IV. Discrete charging states 
accompanied by oscillator excitations. V. Direct mapping to 
CPB. VI. "Pure" charging states. 



of the charging states where the degeneracy in the cross- 
ing points is lifted by the tunneling between the ground 
states in each potential well, Ao being the tunneling am- 
plitude. 



IV. QPS-TRANSISTOR 

The physics of the QPS-transistor is determined by 
the same energy scales as those of the QPS-box: induc- 
tive El, charging Ec, and phase-slip Eg energies. The 
relations between the scales, as outlined in the end of 
the previous Section, determine the qualitative features 
of the QPS-transistor. In this Section, we thus concen- 



trate on the features of QPS-transistor that are different 
from QPS-box or just do not exist there. As mentioned, 
the levels of QPS-transistor are sensitive to flux: the fea- 
ture absent in QPS-box. A symmetric QPS-transistor 
{Es.i — Es,2, El^i — El 2) exhibits a peculiar separation 
of quantum variables that leads to double degeneracy of 
the levels at <£> = ir. We present the numerical results at 
intermediate impedance. 



A. Second-order corrections 

As noted in Section HH the phase-slip corrections to 
the energy levels of the QPS-transistor are of the second 
order in the phase-slip amplitudes. The non-perturbed 
states are labeled with the phonon number n and the 
winding number p. If we restrict the external phase $ 
to the interval [— tt, it] the ground state corresponds to 
|0,0). The correction to the ground state energy reads 

(2) = v |(0,0|gs|n,±l)| 2 
3 ^ E - E ±1 - faon 

E p being inductive energies. To evaluate this expression, 
we represent it in the form of an integral over an auxiliary 
variable t, 



E 



dt 



11 



ei E - El -hu: n)t . |( ,0|#skl)| : 



= (B -S_i-fia;oJV)t 



1(0,01^,-1) 



After this, the sum over n can be taken and the integral 
can be evaluated. The result reads 



4 2) 



E S,l 



„ 2 „, 2 I 2 



^ 2 „, 2 ; 2 
_7r 7 »2 . 



■k 2 {E l>1 +E l>2 ) 

E S,2 

ir 2 {E L ,i+E L a) 

1Es,iEs,2 , . , 
— — - cos 7T- I e 



— 1F1 (a + 7r 2 7 2 ; 1 + a + ir 2 j 2 ; l\ir 2 ^ 2 ) + (a + -> a ] 
a + 



— 1F1 (a+ + tt 2 7 2 ; 1 + a+7r 2 7 2 ; ? 2 tt 2 7 2 ) + (a+ a_) 



9 9 / 1\ 1% 



" + 



1F1 (a+vr 2 7 2 ; 1 + a + ?r 2 7 2 ; -I^tt 2 -/ 2 ) + (a+ ->■ a_) 



(20) 



We have also introduced the notations l\, I2 for the ra- 
tio of the inductive energies and a± for the dimensionless 
energy differences: 

El,i , El.2 



h 



a± = 



, h = 

El,1 + El,2 El,1 + El,2 



n 2 (E L .i+E L: 



The correction Eq. [20] naturally separates into two 
parts: the "classical" one, that is contributed by the 



squares of the phase-slip amplitudes in two wires, and 
the interference one, that bears charge dependence. Its 
rather complicated expression has simple asymptotes at 
large and small 7. 



In the limit of large impedance (7 <§; 1), the correction 
is determined by the contribution of the virtual state with 



11 



n = and reads 



E S,2 



E S,1 



■2E S ,lE S ,2 cos (Trf) 



E 



L.2) 



(21) 



The " classical" and interference part is of the same order 
of magnitude. 

In the opposite limit (7 3> 1), the "classical" part re- 
mains of the same order of magnitude while the charge- 
sensitive interference correction is exponentially sup- 
pressed 



E S,1 



7T 2 (£ 



L,l 
E S,2 



E 



L.2 



1 



1 



TT 2 (E L ^ + E L>2 ) 
2 E S,l E S,2 



a + + If a_ + If 
1 1 



a+ + Z 2 V a _ + ll) 



^ 2 i£L.., -rEL,2) C0S vl)^ hl 2 il - |<T> /7T) 



2tt 



e -X_[l + 2i 1 ;2((l-|*|/7r)(ln(e(l-*/7r))-l)]_ 



(22) 



Interestingly, the exponent of this suppression varies with 
$ changing between tt 2 ^ 2 /2 and 7r 2 7 2 (l — 2lil 2 )/2. 

The correction diverges near the degeneracy points 
<I> = ±7r. This divergence is trivial, indicating the mixing 
of two crossing states |0, 0) and |0, ±1) by the phase-slip 
amplitude. The correction to the ground state near the 
divergence point reads 
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FIG. 15. QPS-transistor: the second-order correction to the 
ground-state energy (in units of E\ = —E 2 S1 / (4(Z?l,i + El,2)) 
and respectively E 2 = -~E s ,iE s ,2/{^{El,i + E Li2 )) ). (a) 
"Classical" part of the correction versus external phase dif- 
ference <E>. From the uppermost curve downwards: 7 <C 1 
(asymptote), 7 = 0.5, 7 = 1, 7 = 2, 7 = 3, 7 > 1 (asymp- 
tote), (b) Charge-sensitive ("interference") part of the cor- 
rection versus "!>. From the uppermost curve down: 7 < 1 
(asymptote), 7 = 4, 7 = 5.65 and 7 = 8. Note the log-scale 
of the plot. 



6E g = -v/e 2 /4+|M| 2 + e/2, 

where e = ir 2 (E LA + E L . 2 )(1 - |$|/tt), and M is the 
matrix element of Hs.i + Hs. 2 between the degenerate 
states, 



M 



-E s ,ie 



i-nqj e 



Es.2 e ~ 



(23) 



This matrix element is exponentially suppressed in the 
limit of small impedance. 

We present the numerical estimations of the second- 
order corrections in Figs. Hnia,b. We assume the sym- 
metric QSP-transistor. The "classical" correction at 
any value of the impedance takes values between the 
two asymptotes, Eq. [5T] and Eq. [2U At each value 
of impedance, the correction diverges upon approaching 
<I> = ±7r. However the divergence controlled by the ma- 
trix element Eq. [521 is exponentially suppressed at large 
7. This is why the curves approach the non-divergent 
asymptote Eq. [22] in the limit of small impedance. 

The charge-sensitive interference contribution is plot- 
ted in Fig. [T5]b in logarithmic scale to show the exponen- 
tial dependence at large 7. Since the exponent depends 
on $, the correction also displays exponential flux sensi- 
tivity. 



B. Large Es 

Let us turn to the QPS-transistor in the opposite 
regime of large Es where the phase-slips produce the 
well-developed potential wells that house the charging 
states and oscillators. As we have seen in the case of 
QPS-box, this situation is achieved for sufficiently large 
E$ for any impedance. 

We concentrate on the QPS-transistor Hamiltonian in 
the charge representation, 



H 



e l,i d 2 E L o ( d 



4 dQ 2 4 \dQ 2 
2E s ,i cos (2ttQi) - 2E S ,2 cos (2nQ 2 ) 
E c (Q 2 -Q 1 ~q/2e) 2 . 



(24) 
(25) 



where the wave function is subject to the periodicity con- 
dition 

*(Qi + i,Q 2 + i) = tf(Qi,Q 2 ). 

The phase-slip terms pinpoint the low-energy wavefunc- 
tions to the minima of the oscillating potential Qi,2 = 
Ai.2- Each minimum corresponds to a charging state 
and gives two series of oscillator levels labeled ni 2. The 
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states are thus labeled as \Ni, N 2 ,ri\, n 2 ) and their ener- 
gies are given by 

E(N h2 ,n 1>2 ) = -\E s ,i\ - \E S . 2 \ + E C {N 2 - JVi - q/2e) 2 
+ hiuirii + hu) 2 n 2 (26) 

with the oscillator frequencies being ftu)\ >2 = 
4tt\/ Es y i, 2 E Ly i^ 2 . 

Interestingly, for a QPS-transistor we get in this limit 
two oscillators instead of a single original one. 

Let us take into account the tunneling between the 
potential wells and do it for the lowest states n\ , n 2 — 
only. The tunneling amplitudes in the directions Q\, Q 2 
are Ai, A 2 e l * respectively, the latter factor incorporates 
the external phase. The expressions for Ai j2 are identical 
to that of the QPS-box, 




FIG. 16. Ground state energy of the QPS-transistor versus 
external phase $(7=1) for a: q/e = and b: q/e — 1. In 
both panels, the phase-slip amplitudes take values: Es/fuoo = 
0, 0.5, 1 and 2 from the uppermost curve to the lowest one. 



A x =8(£| 1 £ L , 1 ) 1 / 4 -exp - 



t V Ela / ' 



A 2 = 8(£ | A2 ,V«. exp ^^j. 
The tunneling part of the Hamiltonian is thus 

i?tunn = - (Al|JVl,JV2)<JVl + 1,^1 

Ni,N 2 

+ A 2 e i9 \N 1 ,N 2 ){N 1 ,N 2 + 1\ + h.c) (27) 

Let us now use the periodicity condition that allows 
us to identify the states with N 2 — Ni — N and disre- 
gard higher oscillator states assuming fiwi, 2 3> Ai, 2 ,i?c- 
With this, the Hamiltonian reduces to 





FIG. 17. Supercurrent in QPS-transistor versus $ for a: 
q/e — and b: q/e = 1. In both panels, the phase-slip am- 
plitude takes values Es/Hu>o = 0,0.5,0.75, 1 and 2. Es = 
corresponds to the straight line. 



power-law suppression in terms of the system parameters, 
I c ~ E 2 j/E c . 

In the opposite limit of Ec 3> Ai )2 the flux sensi- 
tivity of the device is determined by the second-order 
corrections to the ground state, that is, ~ A 2 /Ec- We 
assume the ground state to correspond to N = and 
— 1 < q/e < 1. The second-order correction then reads 



H = V] (-(Ai + A 2 e l *)\N)(N + + h.c.)+E c {N-q/2ef 



N 

(28) 

This is the Hamiltonian of the CPT where the Josephson 
energies Ej,i,Ej, 2 are replaced with 2Ai,2. This reduc- 
tion is less trivial than that in the case of QPS-box since 
QPS-transistor has in principle more degrees of freedom. 

In the limit of Eq <C A12, the QPS-transistor is 
thus equivalent to a double Josephson junction with the 
ground-state energy — |Ai + A 2 e l *| and the critical cur- 
rent I c ~ (2e/fi)min(Ai,A 2 ). 

This estimation should be compared with the maxi- 
mum current through the wire in the ground state under 
the conditions of phase bias. Since the winding number 
in the ground state adjusts itself to the phase, the effec- 
tive phase difference never exceeds tt and this current is 
given by I c {2e/h)'KE Ltl E Ll2 / \E L ,x + E L<2 ) (this is much 
smaller than the critical current in a current-biased wire 
where large phase differences can be built) . Since the tun- 
neling amplitudes Ai. 2 are exponentially suppressed (Eq. 
??) the critical current in the wire is also suppressed ex- 
ponentially in system parameters, ln/c ~ — yJ(E s / El). 
Notably, this suppression is more efficient than the in- 
tercepting the wire with tunnel barriers, that leads to a 



(2) 2|A 1 + A 2 e^| 2 
E c (l-(q/2ey 



(29) 



This diverges at q/e — ±1 where the energies of two 
charging states cross. Near the avoided crossing, the 
ground-state energy is given by 



E g = E c /4- Je 2 + \At + A 2 e 



i$|2 



(30) 



where e = Ec(l ~ q/e) -C Ec (for the crossing at q/e = 
!)• 



C. Degeneracies in a symmetric QPS-transistor 

The odd integer values of q/e correspond to the double 
degeneracy of the pure charging states, while the pure 
flux states are degenerate at half-integer values of the 
external flux $ = ±7r. We have seen that in general 
the degeneracies of this type are lifted: Eq. [5T] and Eq. 
1501 describe avoided crossing of flux and charge states, 
respectively. However, the non-diagonal matrix element 
in both expressions vanishes if the degeneracies were to 
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occur simultaneously, that is, at $ = ±7r and q/e = ±1, 
and the QPS-transistor is completely symmetric (Es,i = 
Es,2,El : i — El,2)- The double degeneracy thus persists 
in this point. 

Since both expressions are perturbative, one could 
think that the degeneracy lifting is just not visible being 
governed by the next-order perturbation terms. However, 
this is not the case: the double degeneracy is preserved by 
a specific symmetry of the QPS-transistor Hamiltonian 
that takes place at odd values of q/e and is not obvious 
from the forms of the Hamiltonian that we have used so 
far. 

To reveal the symmetry, let us rewrite the Hamilto- 
nian of the QPS-transistor in the basis of unperturbed 
wavcfunctions \n,p), n being the number of quanta of 
the original oscillator and p being the winding number 
across the wires. Actually, this is the representation we 
have made use of in numerical simulations. It reads 

H = ^2{E P + haj n)\n,p)(n,p\ 

n , p 

+ J2( P ^n'\n,p)(n',p+l\ + h.c), (31) 



E p = (2Trp+$) 2 E LA E L , 2 /4:(E LA +E L ,2) being the induc- 
tive energies of the states with different winding numbers. 
The matrix elements of the phase-slip terms read 

Pn,n> = -E S ,x^ q,e I dtfrMMn'tt + ^h) 



E. 



S,2 



#^„(<MM0-27rZ 2 ), (32) 



ipn(<fi) being the oscillator wave functions in the phase 
representation. It is evident from this expression that for 
a symmetric QPS-transistor the matrix elements obey a 
specific selection rule: they are non-zero only between 
the oscillator states of different parity, n + n' should be 
odd. Therefore, the whole Hilbert space is separated in 
two disconnected blocks. In one block("+"), n is even for 
even p and odd for odd p, while in another block ("-") 
the situation is reversed: n is odd for even p and even for 
odd p. 

The energies in each block are 47r-periodic functions 
of the external phase. Since 2ir shift of the phase shifts 
p by 1 and therefore switches between even and odd p, 
these energies must satisfy E + (<fi) = E~(2tt + <f>). This 
implies the double degeneracy at all half-integer values 
of the external flux provided q/e is odd-integer. 

In fact, the even-integer values of q/e are also specific 
for the symmetric QPS-transistor. As one can see from 
Eq. [321 a t these values the matrix elements are only non- 
zero between n of the same parity. The Hilbert space sep- 
arates into two blocks comprising respectively odd and 
even n. The energies remain 2tt periodic, so this sepa- 
ration does not imply extra degeneracies at half-integer 
values of the external flux. However, the "random" de- 
generacies due to the crossing of the levels of different 
blocks at some non-specific values of <fi do occur. 
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FIG. 18. QPS-transistor in the limit of large E s . The 
positions of the potential minima in extended charge space 
{Qi-iQi)- The dotted lines connect equivalent minima, rep- 
resenting the set of points which are the same in this space. 
Dark dots give unique minima. Tunneling amplitudes through 
potential barriers (Hamiltonian Eq. 1271) are indicated by 
Greek letters. 




FIG. 19. Flux dependence of the lowest three energy levels for 
a "symmetric" QPS-transistor (El,i = El,2\Es,i = Es,2), 
for 7 = 1, Es/fkjJo = 0.5 and a: q/e = and b: q/e = 1. Note 
the degeneracies at half-integer flux for q/e = 1. 



Although these peculiar separations hold for a sym- 
metric QPS-transistor only, they may become important 
in the course of quantum manipulation of the states of 
QPS-transistor. It is not excluded that the values of Es 
in a QPS-transistor may be tuned by some extra gates, 
this would be a fully symmetric QPS-transistor practi- 
cal. This is why we have discussed it and give numerical 
results for the symmetric QPS-transistor configuration. 



D. Flux sensitivity 

We illustrate the flux sensitivity of QPS-transistor at 
different values of phase-slip amplitude Eg = Eg t i = 
Eg,2- We restrict these examples to the case of interme- 
diate impedance 7 = 1. 

In Fig. \W[ we give the dependence of the ground state 
energies on the external phase $ for weak, moderate and 
large phase-slip amplitudes. For vanishing Eg, the de- 
pendencies are parabolic. For q/e = (Fig. [THla), we 
see the rounding of parabolas at $ = ±7r so the curves 
approach the cos shape upon increasing Es- Also, the 
flux sensitivity defined as the energy difference between 
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FIG. 21. The Josephson circuit dual to QPS-box. 



FIG. 20. Flux sensitivity versus phase-slip amplitude for two 
values of charge (7 = 1). The crosses: q/e = 0. X-symbols: 
5 = 1. 



relation: [Q, </>] = —i. With this, we can readily establish 
the Hamiltonians of the dual devices. The Hamiltonian 
dual to that of QPS-box Eq. [3] reads 



<J> = 7r and $ = 0, quickly decreases upon increasing 
Es- For q/e = 1, the curves remain parabola-like, this is 
due to special degeneracy described in the previous sub- 
section. Nevertheless the flux sensitivity decreases upon 
increasing Es, although not as fast as in q/e = case. 

The energies of the three states with lowest energy are 
plotted in FigH91 at Es = 0.5hw - At q/e = 0, this ex- 
emplifies the degeneracy lifting and avoided level-crossing 
at half-integer flux. In contrast to this, we see unavoided 
crossing of the first and second excited state at $ sa ±2.3, 
this indicates opposite parity of the states. At q/e = 1, 
we see the double degeneracy at $ = ±n. 

Fig. [2TJ1 gives the flux sensitivity versus Es for q/e = 
and q/e = 1. They coincide at Es = indicating no 
charge sensitivity. Most interesting region is that of mod- 
erate Es where the flux sensitivity is charge sensitive but 
still substantial. In agreement with the considerations in 
Subsection B, the flux sensitivity drops exponentially at 
large Es. Although this takes place for both values of the 
external charge, the flux sensitivity drops much slower for 
q/e = 1 and eventually exceeds by far the flux sensitivity 
at q/e — 0. Indeed, the same considerations show that 
the phase-dependent energy is ~ A for q/e — 1 and ~ A 2 
for q/e = 0. 



V. DUAL DEVICES 

As discussed in2£, there is a duality between phase-slip 
junctions and Josephson junctions. Each device contain- 
ing phase-slip junctions has an analogue where phase- 
slips are replaced with Cooper-pair tunneling events, 
phase is replaced with charge. Upon this dual transfor- 
mation, the devices have identical quantum dynamics. 
In this Section, we will find and shortly discuss the dual 
analogues of the QPS-box and QPS-transistor. 

To avoid any misunderstanding: the analogies be- 
tween CPB and QPS-box, and between CPT and QPS- 
transistor that we have thoroughly discussed in this ar- 
ticle are not related to the duality under consideration. 
Dual devices are rather different. 

The duality is based on the canonical transformation 
(Q,4>) — >• (— 0/27T, 2ttQ) that preserves the commutation 



# QPS-bo* = ^ (4> ~ $) 2 + E' C Q 2 - Ej cos I (33) 

It is a Josephson Hamiltonian with the parameters re- 
lated to the initial QPS values as follows: 

Ej -> 2E S ; E' L -> -f ; E' c -> ir 2 E L - $ -> nq/2e 

In electrical terms, this is a Josephson junction in series 
with an LC-oscillator. The flux sensitivity of this device 
is obtained from the charge sensitivity of the QPS-box, 
E($) = £ QP s-box(<7 = e$/7r). 

For the QPS-transistor the same transformation leads 
to the following dual Hamiltonian: 

^QPS-transistor = ^i(<fe - 4>l ~ $) 2 

+ ^ i g 2 --Bj,iCos(^i) 
+ (!)-> (2). 

where the constraint Qi + Q2 = q/2e + P is imposed, p 
being an integer number. Physically, this is an integer 
number of of Cooper pairs accumulated in the nodes en- 
circled in Fig. [221 We note that in the limit of vanishing 
Ej this combination of nodes is isolated from the leads 
and thus can sustain integer charge. To keep consistent 
notations, the gate capacitance in Fig. |2"21 should be as- 
sumed vanishingly small, C g -C C1.2, whereas C g V g = q. 

The regime of vanishing phase-slips in original devices 
correspond to the regime of vanishing Josephson cou- 
plings for dual ones. In this regime, the devices are 
LC-oscillators. The dual of QPS-transistor possesses an 
extra degree of freedom: the number of integer charges 
accumulated in the isolated island. The regime of well- 
developed Coulomb blockade in the QPS devices corre- 
sponds to classical limit of the dual Josephson circuits, 
where the winding numbers in the superconducting in- 
ductors are dual analogues of discrete charges. 

VI. CONCLUSIONS 

In this paper we have proposed and discussed two 
superconducting devices made of superconducting wires 
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FIG. 22. The Josephson circuit dual to QPS-transistor. 



subject to coherent quantum phase-slips: QPS-box and 
QPS-transistor. Our main goal was to demonstrate the 
charge sensitivity, this would be unambiguous experi- 
mental signature of Coulomb-blockade behavior. The ex- 
perimental realization of our proposed devices is achiev- 
able with the state-or-the-art technology. 

We have shown that the charge sensitivity appears al- 
ready for small phase-slip amplitudes as a perturbative 
correction to the ground-state energy. This correction 
is of the first order in case of QPS-box and of the sec- 
ond order in case of QPS-transistor. In both cases, the 
charge-sensitive part of the perturbative correction is ex- 
ponentially suppressed in the limit of low impedance. 
In contrast to the QPS-box, the QPS-transistor exhibits 
both flux and charge sensitivity that makes it potentially 
useful for measurements. 

However, if the phase-slip amplitude becomes of the 
order of either charging energy (large impedance regime) 
or inductive energy (small impedance regime), both de- 
vices show discrete charging states that follow a common 
Coulomb blockade pattern of " crossing parabolas" . The 
crossover to Coulomb blockade occurs differently in the 
limits of large and small impedance. We have revealed 
and investigated both analytically and numerically six 
distinct parameter regions. For QPS-transistor, we have 
analyzed the flux sensitivity (superconducting current) as 
well as combined flux-charge sensitivity . For a symmet- 
ric QPS-transistor we have found a variable separation 
at specific values of induced charge that leads to double 
degeneracy of the states at q/e = 1 and half- integer exter- 
nal flux. We have calculated the superconducting current 
through QPS-transistor to show the non-triviality of the 
device. We briefly discuss the Josephson-based devices 
that are dual to QPS-box and QPS-transistor. 
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Appendix A: Microscopic foundation 

In this Appendix, we sketch the microscopic reasoning 
that justifies the Hamiltonians {3}, ([5]). Most of this rea- 
soning can not be regarded as the resultive part of the 
article. Rather, this Appendix mainly summarizes the 
results of Refs . 15 i 23 i 24 i 26 perhaps with more comprehen- 
sive notations. We include this summary for the sake 
of completeness, and to respond to the requests of our 
colleagues. 

The reasoning proceeds in two steps. At the first step, 
one recognizes how the phase-slip amplitude emerges 
from instantons of an effective field theory describing the 
quantum fluctuations of the superconducting order pa- 
rameter. The outcome is a Id sine-Gordon-type model 
of the superconducting wire where the phase-slip am- 
plitude per unit length enters as a parameter. At the 
second step, one replaces this Id model with a zero- 
dimensional one. This is justified if the internal degrees 
of freedom of the wire can be efficiently disregarded or 
incorporated into the parameters of the zero-dimensional 
model. We discuss the requirements for this. The out- 
come is a Hamiltonian combining inductive energy of the 
wire and total phase-slip amplitude. 

One can start with a Hamiltonian of electrons 
in the metal wire subject to Coulomb and phonon- 
mediated interaction. Hubbard-Statonovitch transform 
in imaginary-time^ represents these two interactions in 
terms of the (quantum) fluctuations of two fields: V(r, r), 
time-dependent voltage, and A(r,r), complex super- 
conducting order parameter. Electro-neutrality condi- 
tion enforces Josephson relation 2eV = h(f> in each 
point and time moment, <fi being the phase of the or- 
der parameter—. After that, one can integrate out the 
electron degrees of freedom and obtain the effective ac- 
tion for A(r,i)i£. If the lateral dimensions of the wire 
are smaller than the superconducting coherence length £, 
it suffices to consider only the realizations of the order 
parameter that are constant along the wire cross-section. 
Those can be parametrized with a Id field A(x,r). The 
action consists of two terms ^ 

„ „ / h\ 2 C f , dr - 2 , . 

s = Ss» P +[y e ) ¥ J dx ^ (x ' r) 

Here, the first term incorporates everything related to 
superconductivity while the second term represents the 
charging energy of the wire, C being the wire self- 
capacitance per unit length. 

The superconducting part of the action is non-local 
in both coordinate (at the scale ~ £) and imaginary 
time (at the scales H/Aq, Ao being the saddle-point value 
for the superconducting gap). The degenerate topolog- 
ically trivial saddle-point solutions correspond to con- 
stant modulus of order parameter and are parametrized 
by the phase, A(x,t) = A e^. This suggests the im- 
portance of Goldstone modes that are long-wave topo- 
logically trivial fluctuations of the phase. The effective 
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" hydrodynamic"— action for these fluctuations reads 

*-®'M<£) + (§*) (A1) 

L being the wire inductance per unit length that is mostly 
contributed by kinetic inductance of the superconducting 
material. This action describes the propagation of the 
electromagnetic waves along the wire with velocity cq = 

i/Vlc. 

The full action also permits topologically non-trivial 
saddle-point configurations, instantons, those correspond 
to phase slips. A single phase-slip configuration has a 
core at a point (xo,to) of a size £,h/Ao, this size cor- 
responding to the non-locality of the action. Far from 
the core, the modulus of the order parameter equals Aq, 
while the phase is given by 

(f>(x, r) = — arctan [ — ^— - ) 

\co(t-t )J 

There is a significant extra action associated with a phase 
slip, Sqpc- This action contains "hydrodynamic" part 
that comes from the long-wavelength action (Eq. IA1|) 
and "core" part. With logarithmic accuracy, the hydro- 
dynamic part is given by 

c nh 2 [C. (x Lc \ 

where -?Q. C . is the lower cutoff (of the order of system 
size) and X u . c . is the upper cutoff (of the order of core 
size). The core part of the action is readily estimated as 

Sc ' QPS = G^M 

where R is the normal-state wire resistance per unit 
length, and a is a dimensionless coefficient which is not 
known so far. The evaluation of the partition function 
of the system in the vicinity of the phase-slip saddle- 
point solution gives the QPS amplitude per unit length, 
Es oc e~ SQPC / h . The total amplitude has been estimated 
inl£, where it has been called "rate" Jqps- We should 
perhaps stress that by no means such jqps gives an esti- 
mate of a transition rate of any event involving quantum- 
phase slips, nor this was indented by the authors of—. 
The transition rate in any case should be proportional 
to 7qp5, the square of the transition amplitude. Un- 
fortunately, the fact that the coefficient a is not known 
makes the theory quite useless for predictions of any con- 
crete values of Eg: The exponential dependence on a 
would result in uncertainties of many orders of magni- 
tude. The core part of the action has been found for a 
ballistic single-channel wirej^i but this hardly helps in 
any realistic situation. 

Further analysis presented i n 15 i 27 concerns the instan- 
ton configurations with multiple phase-slips. The prob- 
lem can be mapped onto a 2d gas of charged particles 



with logarithmic interaction. For our purposes, it is con- 
venient to rewrite it as an equivalent sine-Gordon model. 
A convenient variable for this is the charge 2eq(x, r) 
passed thought the point x of the wire (the charge density 
accumulated in the wire is then —2eq' while the current 
I = 2eq. The action in terms of this variable takes the 
following form 

-E s cos(2Trq(x 1 t))) (A2) 

with two first terms giving the hydrodynamic part of the 
action, that is, (|A1|) . The quantity Es corresponds to 
"bare" fugacity of the phase-slips, that is, contains only 
the "core" contribution of the action, Es oc e~ Sc > QPS / h . 
The uncertainty in defining the core is incorporated into 
upper cutoff of the sine-Gordon model. 

From a formal point of view, it seems natural to con- 
sider an infinite wire. In this case, it is very well known 
that the system exhibits a quantum Schmid transition 
(sometimes termed BKT, from analogy with associated 
2d classical system) at arbitrary Es and at a critical value 
of the wave impedance 

IT _ h 
\I~C~ 8e?' 

In the context of the wires, it is sometimes called 
the superconducting-insulating transition. Perhaps one 
should be more cautious about this, since the essence of 
the transition is the distinct behavior of the renormal- 
ized inductance L of the wire as function of the wire 
length I, L(l). The renormalization by the phase-slip al- 
ways enhances the inductance, the enhancement changes 
from power-law L oc l& with /3 < 1 to exponential depen- 
dence L oc e l l lc at the point of the transition. Guided 
by personal taste, one may call the wire insulating (zero 
inductance per unit length in the limit / — > oo) or super- 
conducting (finite inductance at any finite I). 

More insight into the wire properties in the insulat- 
ing regime can be obtained when considering the saddle- 
point solutions of the sine-Gordon model. The trivial so- 
lution corresponds to q = so that no charge can pass the 
wire. Besides, there are sine-Gordon kinks. In this case, 
they correspond to the charged soliton-like excitations. 
The minimum energy of the soliton, E so i ~ eyEs/C , is 
the gap for charged excitation. The soliton is spread over 
the length i sol = e/^E s C. The inverse inductance, that 
is, the supercurrent in the wire, is due to tunneling of 
the solitons over a potential barrier of the length I and 
the height E so \. Therefore, the wire in the insulating 
regime is a kind of a semiconductor, a uniform barrier 
for charges to tunnel through. The ground state is ho- 
mogeneous, one cannot say that the wire is cut into small 
pieces separated by tunnel barriers. 

It is not clear at the moment to which extent the Id 
model can be useful to describe experimental results even 
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at qualitative level, this being a subject of ongoing re- 
search and debate. The realistic wires are not only fi- 
nite: they are rather short. The inverse time of electric- 
ity propagation across the wire, cq/1, can easily exceed 
the typical superconducting energy scale Aq/H. In this 
case, the upper cutoff is smaller than the lower cutoff, 
this renders hydrodynamic corrections irrelevant. Even 
if this is not the case, it has been showr* 2 ^ that for a finite 
wire embedded into an external circuit the Schmid tran- 
sition should be governed by the external impedance of 
the circuit rather than the wave impedance of the wire. 
There were recent attempts^ to modify the traditional 
rcnormalization schemes to take into account finite I and 
possible normal excitations not captured by the hydro- 
dynamic action. More experiments and more detailed 
comparison of experiment and theory are required to re- 
solve the issue. 

A way to avoid these complicated issues is to think 
of a phase-slip superconducting wire in terms of a zero- 
dimensional model, that is, in terms of a circuit-theory 
element. For a closed wire, such model has been proposed 
inSl. The model can be derived from the sine-Gordon 
model (|A2|) by setting q = const (x). The result is most 
convenient to present in the Hamiltonian form 

H = El<\? — E s cos(2irq) 



El being inductive energy of the wire, E$ = Egl. The 
operators of phase drop across the wire, <j> and charge 
passed q satisfy canonical commutation relation [q, cj)] = 
i. This zero-dimensional model has been extended in 26 
to the case of the wire embedded into an arbitrary exter- 
nal circuit and its exact duality with Josephson-j unction 
circuits has been demonstrated. 

The zero-dimensional model obviously does not de- 
scribe the internal excitations of the wire such as standing 
electromagnetic waves in the absence of the phase slips 
or charge solitons arising from sine-Gordon model in the 
limit of sufficiently developed phase-slips, although the 
renormalization of the parameters El and E$ by these 
excitations can be taken into account on phenomenolog- 
ical level, that is, just by taking the measurable val- 
ues of these parameters rather than the "bare" ones. 
One can argue that these excitations either have energies 
exceeding the energy scales of zero-dimensional model 
Es,El ^> %cq/1,E so \ or just do not fit into the wire 

(Aral > 0- 

A convenient formal way to assure irrelevance of in- 
ternal excitations is to take the limit of vanishing self- 
capacitance of the wire, C — > 0. In this limit, all of the 
above requirements are fulfilled, since the energy scale 
of the standing waves, soliton energy and length become 
infinite. The charge cannot be accumulated in the wire, 
q'(x) = 0, and the model is essentially zero-dimensional. 
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